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We use the exact correlation function of the density of energy levels in the magnetic field to 
evaluate persistent currents in mesoscopic metals. We also analyze the perturbation theory limit of 
the correlation function vis-a-vis the perturbation theory limit of the orbital response. 

I. INTRODUCTION 

Orbital magnetism and persistent currents in mesoscopic metal rings and grains have received much attention in the 
past decade. Theoretically, mesoscopic orbital magnetism is rather well understood for temperature T> A, where 
A is the mean energy level spacing, and for magnetic fields such that the flux through the sample i^] <j> < 0o = hc/e 
- the flux quantum |lj (<fio — 2n/e in units where h = c = 1). To evaluate the thermodynamics quantities at these 
temperatures, it is sufficient to characterize the energy spectrum in terms of the level density p. In particular, the 
orbital magnetic dependence of the level density correlation function, 

R(cj,T H - 1 ) = (5p(e 1 ,H)5 P (e 2 ,H)) (I) 

where dp [e, H) — p (e, H) — (p) and the brackets stand for disorder averaging, describes both the mean mesoscopic 
response and the fluctuations |lj . 

Furthermore, in this limit it is sufficient to use the perturbative limit of the correlation function. In the two- 
Cooperon approximation, it is given by 

s 2 1 

R pert (w,t» ) = . (2) 

u) = (e2 — £i), 7 is the level broadening (which will be neglected below relative to T, as it should for metals), and s 
and [3 are integers whose values depend on the symmetry class. For Gaussian Orthogonal Ensemble (GOE) (3 = 1, 
while for Gaussian Symplectic Ensemble (GSE) j3 — 4 due to the fact that 3 out of 4 spin channels are eliminated by 
spin-orbit scattering. Neglecting Zeeman splitting, s — 2 in either case due to double spin degeneracy and Jo) A = s. 
th^ 1 gives the orbital magnetic dependence. For a single period for a ring with the flux cf>, it is given by pi 

T H ^ = {2,?E c (^j (3) 

where E c — HrrJ^f 3> T and T<nj f — L 2 / D is the time of diffusion around the ring of circumference L (time of diffusion 
to sample boundaries for simply connected geometries, such as disk D being is the diffusion coefficient. 

In the evaluation of the mean response one must remember that the electrons in a closed system form a canonical 
ensemble |2| . Ingeniously, Imry || found a simple relationship between the free energy of a canonical ensemble and 
the free energy and particle- number fluctuation of an equivalent grand canonical ensemble (see Appendix A) . It was 
shown Q], however, that such a relationship holds only if the particle-number fluctuation is large, which is the case 
when T ^> A also. A perturbative evaluation JjJ using Imry's formalism yields a persistent current oc T — 1 . The 
validity of this result is obviously limited to temperatures T > A. 

Evaluation of the orbital magnetic response for T < A is a challenging problem. Unlike spin magnetism, which can 
be analyzed in terms of a few-level problem ||, all states of the Fermi sea contribute, in principle, and there are large 
cancellations between the diamagnetic and paramagnetic components. As a first step, however, one can use a "mixed 
approach" wherein Imry's formalism is combined with the exact form of the level-density correlation function. Such 
an approach should, hopefully, reveal a trend in the persistent current as T becomes smaller than A. 
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1 For a narrow ring, the response is <j>o/2 periodic so there is no need to consider larger fluxes. 
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Two earlier attempts of this nature were undertaken in Refs. M and M. While both correctly argued that the 
apparent singularity for T — > found in the perturbative approximation disappears for GOE as the persistent current 
saturates to a finite value as T crosses A, they bore serious deficiencies. The former did not find a closed form of 
the correlation function and considered only T = (while Imry's formalism, as was pointed out above, is rigorously 
applicable when T> A). The latter considered the entire range of temperatures but the magnetic field dependence 
of the correlation function was crudely mimicked by a substitution 7 — * 7 + th^ 1 (as in perturbative expression (|^)) 
in the exact correlation function obtained in the assumption of presence of level broadening 7 alone. Later, however, 
the " mixed approach" was fully implemented for spin magnetism [Q where the magnetic field dependence of the level 
density correlation function is rather simple. 

It turns out, however, that the exact correlation function in the presence of a time-reversal symmetry breaking field, 
such as magnetic field, had been calculated over a decade earlier by Mehta and Pandey || for both GOE and GSE 
using random matrix theory (RMT) (their result was subsequently confirmed using the supersymmetric non-linear 
a- model in Ref. ||). In this work we combine the exact correlation function obtained in Refs. || and || with Imry's 
formalism to find the magnetic field dependence of the persistent current for finite temperatures. 



II. LEVEL CORRELATION FUNCTION 

In units of (p) 2 , the dependence of the correlation function of the level density on the time reversal symmetry 
breaking parameter A in transition from GOE to GUE (Gaussian Unitary Ensemble) and GSE to GUE is given by 

§: 

^OC pi T,/ 

Rorth (x, A) = / dk sin (kx) sin (k 1 x) e~ 2xV ^ (4) 

Ji Jo k 

-1 pOO pi 1 

R sym p(x,X)= 7 dk dk'— sin(kx)sin{k'x)e- 2x27r2 ( k2 ~ k ' 2 ) (5) 
4 Ji Jo * 

respectively, where x = ^l^i — £2 1 ■ For a ring enclosing the flux <f> ||], 

A 2 = (47rAr H ) _1 (6) 

where th^ 1 is given by eq. (j^). 

It follows from eqs. ([|) and (^) that 

d f sin [x) \ /7T 



^(^0)=-^(^j(|-Si(^)) (?) 
„ , 1 ( , x d ( sin(x)\ , A 

Rsymp (X, 0) = - ( TT5 (x) + — ( — ^ ) Si (x) ) (8) 



The main difference between GOE and GSE is the (S-function term in R sym p {%, 0). More precisely, it is a A — > limit 



The latter will prove to be responsible for the singular temperature dependence of the persistent current in GSE in 
the "mixed approach." 

Consider now the full correlation function of the level density Q , 

sin 2 ( x} 

n orth {x, A) = nS (s) ±-L + R orth (x, A) (10) 

x z 

1 / sin 2 (x) 

and notice that for either ensemble R(x,X) — > when A — * oo. It is clear then, using eqs. (0) and (§), that GOE, 
GSE, and GUE all exhibit a property of "level rigidity," as manifested by 

J K (x, 0) dx = J K (x, oo)dx = (12) 

(in fact, it is likely that J 1Z (x, A) dx = for any A). 



Tisyrap (x, A) = - 7T(5 (x) g ) + R symp (x, A) (11) 
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III. PERSISTENT CURRENTS 



As was explained in Refs. (^] and [[I), the dominant contribution to the mean orbital magnetic response of isolated 
metals is given by the part of the magnetic free energy which differentiates a canonical ensemble from an equivalent 
grand canonical ensemble (whose mean part icle number coincides with the fixed particle number of a canonical 
ensemble). For T A, it is given by eq. ( |A9| ) of Appendix A, 

SFfI = 2jp) dH J J /( ei '^)/( e 2'^)- R ( e i' £ 2,TH _1 )deid£ 2 

1 ] h I dio -Rlu, th- 1 ) (13) 



2(p) "J exp(f)-l 

where integration on {e\ + £2) was performed. Once SFh is known, the persistent current and the magnetic moment 
are obtained, respectively via the following thermodynamic identities 

We will use the expressions for the correlation function given in the previous Section to evaluate eqs. ( |l3| ) and (|l4|). 



A. Perturbative limit 



Here we address how the perturbation theory limit, whose leading term is given by eq. (^|), is recovered from eqs. 
(||) and (|^) (as was already mentioned, in what follows we will neglect 7). This limit is important for the evaluation 
of the orbital magnetic response when T> A. The linear-response part, in particular, is defined by the term cx 1 
in the free energy. Notice, however, that the expansion of eq. (0) in orders of t^ 1 yields no linear contribution. The 
contribution linear in in the level correlation function comes, in fact, from the next order term in expansion of 
eqs. (0) and (^) and corresponds to the three-Cooperon diagrams, while it is the two-Cooperon diagrams that yield 
eq. (0). On the other hand, we will see that in the expression for the free energy such an expansion of the correlation 
function is not well defined and a careful analysis shows that the two-Cooperon diagrams give a leading contribution 
to the linear response, while the contribution of the three-Cooperon diagrams is smaller in the parameter A/T. 

Eq. (0) can be written as 



Rorth (x, A) « Im • 



D (ix-4\ 2 Ti 2 )k 



c-4X 2 



-dk > Im 



1 



dk'k'e( tx+4X ^ 2 ) k ' 



(15) 



- 4A 2 7T 2 



(ix - 4\ 2 ir 2 Y 



Jx+4\ 2 ir 2 



1 



1 



4A 2 7T 2 



(ix + 4A 2 7r 2 



c.c. 



where we used the fact that for the large values of \ix — 4A 2 7r 2 | the main contribution will come from k,k' ~ 1, so 
that k 2 w 2k — 1, k' 2 w 2k' — 1 and used integration by parts. Omitting oscillating terms in eq. (|l5|), we find 



R or th (x, A) w ^ Re • 



Recalling that (p) A = s, we finally obtain 

s 2 

Rorth (w, t h ) w — — ■ Re 
2pir z 

A similar calculation yields 



1 



Rsymp {^iTh 



Re 



{-ix + 4A 2 7r 2 ) 


2 Hx + 4A 2 7r 2 ) 3 




2 


[ (-iuj + TH^ 1 


) 2 {—iu + th^ 1 ) 3 


I 1 


2 

H 



2(3ir 2 l(_^ + Tff -i) 2 (-iu + TH- 1 ) 



(3 = 1 



(3 = 4 



(16) 



(17) 



(18) 



The first term in eqs. ( |l7| ) and (18 coincides with eq. (||) and the second term represents the three-Cooperon diagrams. 
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Substituting eq. (||) in eq. jl^), we find 



6F H = ~ W) d H J duR pert {u,r H -) cxp( . ) _ 1 (19) 

Expansion of R per t in terms of th^ 1 would lead to divergence of the integral at zero. Consequently, one needs to treat 
carefully the second order pole at -irn . Alternatively, one can change the order of J and Re and the pole will be 
of the first order. Using the latter procedure, one can close the contour in the plane that does not contain the pole 
and evaluate the integral in terms of the sum over Matsubara frequencies, u) m — 2^imT . Since each summand is oc m 
and rjy -1 is finite at this stage of the evaluation, the m = term yields zero and the expansion in orders of th _1 is 
possible in the m > 1 Matsubara sum. As a result, one finds the formulae derived in Ref. [jjj. The second term in 
eqs. ([l7] ) and (|l^) can be treated similarly and, as was mentioned above, gives a correction of the order of A/T. 

For completeness, a detailed perturbative evaluation of the persistent current is given in Appendix B. Its main 
conclusion is that the persistent current is given by 

s 2 8,(2)eE c <j> 

Ipert {9} ~ rp/ \ 1") 9 < 9c (20) 



2/3tt 2 T ( P ) 0o ' 

s 2 e tf) Q 
2(3ir 2 2tt (p) <j> 



Ipert (</>) « * ■ o C , A 4>c « 9 « 0o (21) 



where 



T 



4nE r 2 



« (22) 



is the magnetic flux scale for GOE to GUE transition. We emphasize the universal character of the decay of the 
persistent current with the flux in eq. ( ^l| ) first noticed in Ref. jj| . 



B. "Mixed approach" 

We now proceed to use the exact level correlation functions (^) and (^|) in eq. (|l|). While the latter is valid for 
T 3> A, in which case the perturbative expression (||) for the correlation function is sufficient, we hope that the results 
obtained in such an approach will be qualitatively indicative of the behavior for T < A. 



1. Orthogonal ensemble 

For the orthogonal ensemble eq. (|l4|) yields, with the help of eq. (^), the following expression for the persistent 
current 



Iorth (9) — T. TTl 



2 d(j> A 

00 pOO pi 1 I 

V m / dk dk'^- (k 2 - k' 2 ) L- 2 - 2 mU^') - e -^l(Ht') e -**<AT*)- l (A a -k") (23) 

m=l Jl Jo k 

where the integral over u was reduced, by means of contour integration, to summation over the Matsubara frequencies 
u> m = 27rmT, m = 1,2, 3... . Performing summation over m, we obtain 

Iorth {9) ~ Y ^ A 

dk f 1 dk 'K. ( fc 2 _ f ^1 ^1 ] g-^CAT*)- 1 ^-*") (24) 

7o ^ V 7 ^sinh 2 (7r 2 g {k-k')) sinh 2 (tt 2 £ (fc + *')) / 

The integral in eq. ( |24| ) cannot be, in general, evaluated analytically and we proceed to analyze some limiting cases 
which can be explicitly handled. 
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a. Perturbative limit, T A While we have already discussed the perturbation theory limit, T S> A , it is 
tutorial to derive the latter directly from eq. (^4|). To this effect, omitting the exponentially small term with (k + k') 
and taking into account the fact that the main contribution comes from k, k' ~ 1, we find 

ww « /* Cm .*-*} e -..^»- ( >->-) 

4 A ocj) J 1 J smlr (7r 2 ^- (fc - A:')) 



^^m^ r*,^^ r du (25) 

J — u 



A 9^ 7 ""sinh 2 (tt 2 ^ 



where we hrst shifted the variables of integration by 1 (since the integrand depends only on k — k') and then introduced 
the new variables u — (k — k') and v = k' and extended the upper limit of integration to oo due to the exponentially 
decaying integrand. Re-scaling the variable of integration, x — (irT/ A) u, we finally arrive at eq. (B2) in Appendix 
B. 

b. Non-perturbative limit, T -C A Expanding the sinh- functions in eq. ( p4[ ) , we find 

<*> » 5TT f * I ' *(P^ e -*" <4 ™ , "' < '"-'" ) < 26 » 

This persistent current is the limiting current in the zero-temperature limit. Notice, that the upper limit of the integral 
on k has been extended to oo (despite the small-parameter expansion of the integrand) , which introduces a correction 

of order max|7r 2 T/A, (i7rr^ 1 /A) 1 ^ 2 | to eq. (|2^). As in the perturbative case, eq. ( f26[) has the following limiting 

behaviors: 

Iorth {(f) » r ^- - 8eS c ^-, « C (27) 

2n o<p tfio 

Iorth M~ o!ao m T-^«^« 00 (28) 



2/3tt 2 2tt (p) 

where now C is found from the condition \kt^ / A ~ f as 



27r 3 £ n 2 

and we used the fact that 



A "°«0 O (29) 



dk / dfc' 



k' 2 1 



(fc 2 - /c' 2 ) 2 

and took explicitly s = 2 in eq. (p7|). Notice that the universal decay of the persistent current with the flux for 
4> c *C 4> ^ 0o coincides with the perturbative result of eq. (pT{): this is obvious both on physical and mathematical 
grounds and comes from fc *~ k' ~ 1. 



2. Symplectic ensemble 

For GSE, we obtain, analogously to GOE, the following expression for the persistent current 

As for GOE, we proceed to investigate the limiting cases of this expression. 

a. Perturbative limit, T A This case is identical, up to coefficient to the orthogonal ensemble since the 
reversed ratio k/k' relative to the latter is inconsequential in view of the fact that the main contribution to the integral 
comes from k ~ k' ~ 1. As a result, we again find that the expression for the persistent current given in Appendix B, 
with the limiting behaviors described by eqs. ( |20| ) and (pl|). 



sinh 2 (tt 2 J (fc + fc')) 



T(Ar„)- 1 (fc 2 -fe' 2 ) 



(30) 
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b. Non-perturbative limit, T <C A As for GOE, expanding the sinh- functions in eq. (|30|), we find 

which is the limiting current in the zero-temperature limit. We first address the zero-field limit. Towards this end, we 
notice that 

which can be easily obtained by differentiation on a and subsequent evaluation of the integral in terms of the Error 
Functions. Consequently, we obtain 



s dr H x 1 / 2A 



hymp (0) « fl £- ■=* — = —>J2^E C , T = (33) 



The latter indicates the existence of the permanent current (and orbital magnetic moment) for GSE, as first conjectured 
by Kravtsov and Zirnbauer |To| . 

The permanent current in a zero flux implies a Curie-type response for non-zero temperature. To obtain the latter, 
we need to be more careful when expanding the sinh- functions in eq. (|3^). Namely, we start by expanding in 
7r 2 ^fc' <C 1 which yields 

where the latter condition on the flux gives the range of applicability of the linear response for GSE in the limit 
T« A. Rescaling the variable of integration and omitting terms small in ratio A/T, we find 



T U\~ 1 9t h 1 a2 (p) r w 3 cosh fa) 2seA v 



From the latter, I sym p {4>) ~ ( s e//3) \/2-k1AE c at ^ ~ (</>o/2) T^/2-Kj (AE C ) - in qualitative agreement with eq. (p3|) . 



IV. DISCUSSION 



The most interesting result is a possibility of a Curie-like orbital response in GSE at temperatures smaller than the 
mean energy- level spacing. In the "mixed approach" considered in this work, it appears as a result of broadening by 
the magnetic flux of a S- function term in the level density correlation function fiof , in accordance with eq. (|]). While 
this broadening was attributed to lifting of the Kramers degeneracy by the magnetic field, the level density 
description of the spectrum is, obviously, insensitive to the differences between even- and odd-electron systems (this 
equally applies to the use of eq. ([l3|) which assumes large variations of the particle number in the equivalent grand 
canonical ensemble). On the other hand, the Kramers theorem distinguishes between even- and odd-electron systems. 
In the odd case, the Curie orbital response would be due to a permanent persistent current (moment) in a state 
"unmatched" by a time-reversed state of a degenerate doublet. The "mixed approach" would then mimic the average 
behavior between the even- and odd-electron systems. The situation is in complete analogy with spin magnetism 
where the "mixed approach" is indicative of the Curie response but one needs to study || few-level systems at 
T <C A to correctly account for odd/even differences. A similar analysis of the orbital response is complicated by 
the need to consider large Fermi sea cancellations between its diamagnetic and paramagnetic components. While 
odd/even difference is not important for GOE, since the Fermi level states are not current-carrying, the few- level 
analysis of the van-Vleck response must be carried out for T« A. We will address these issues in a future work. 
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APPENDIX A: RESPONSE OF THE CANONICAL ENSEMBLE 



In this Appendix we derive the relationship between the free energies of the canonical ensmble and an equivalent 
grand canonical ensemble in the limit T> A. Denoting du the change of the appropriate quantities with the magnetic 
field, we find 



d H F = d H Sl + Nd H p 



(Al) 



where it is assumed that the chemical potential of a system may depend, generally speaking, on the magnetic field. 
Using the idea that the magnetic field dependence can be transferred from the energy level dependence to the density 
of states 12 , 



we find 



where 



= -T J dep (e, H) In ^1 + exp 



-T J dEd H p(s,H)bx(l + exp 



p — e 
T 



N = dep (e, H)[l + exp 



T 



e — p 



(d H p) N 



T 



(A2) 



(A3) 



(A4) 



is the fixed number of particles in the canonical ensemble and the mean number of particles in the equivalent grand 
canonical ensemble. Combining eqs. QA1[) - (A4), we obtain 



(A5) 



d H F = ~T J ded HP (e,H)\n ^1 + exp^^^ 



where the (dnp) N term, containing the dependence of the chemical potential on the field, canceled out. Notice that 
this is in agreement with the small corrections theorem wherein the corrections to F and f2 are the same, but in the 
latter we do not vary p. The derivation that follows is insensitive to the dependence of p on the magnetic field but 
only to the fact that it varies from system to system B. 

Denoting (p) the average value of p over all systems (disorder-average) in zero field, we proceed to expand in powers 
of H - (p) 



d H F « -T J ded H p (e, H) In ^1 + cxp 



</<) 



T 



-(p-(p))J ded H (p(e,H)-(p))f(e,(p)) (A6) 

where f(e,(p}) is the Fermi distribution function and the first term has the meaning of djjQ ({p))- Taking into 
account the following expression for the number of particles, 



N = J de(p)f(e,(p)) 
= J dep (e,H)f(e,p) 

de( P )f(e,(p))+ I de(p(e,H)-(p))f(e,(p))-(p-(p)) I de (p) f (e,(p)) 



(A7) 



and the fact that the last integral is — (p) , wc obtain, comparing the first and the last lines of the above equation, 



( M -( M )) Rj -(p)- 1 / de(p(e,H)-(p))f(e,(p)) 



(A8) 



Substituting the latter in eq. (^9) and performing disorder averaging, we find 
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SF H = d H F - d H Q. ({ft)) 



/ j dsids 2 (5p(e 1 ,H)d H 8p(e 2 ,H))f(si,{p))f(e 2 ,(p}) 



1 



2(P) 
1 

W) 



Oh 



deide 2 (Sp (e%,H) Sp (e 2 , H)) f (s 1} (p)) f (e 2 , (p)) 



(A9) 



d H (SN 2 ) 



where ((57V 2 ) is the mean particle number fluctuation in the equivalent grand canonical ensemble |l3| described by 

n((p)). 

The transformation in eq. ( |A9| ) can be explained as follows. Notice that (5p (si, Hi) dn 2 Sp (^2, H 2 )) can be written 
as 8h 2 (Sp (eij-Hi) Sp (e 2 , H 2 )). However, the diffuson contribution to (Sp (e\,H\) Sp (e 2 , H 2 )) depends on (Hi — H 2 ) 2 
and the Cooperon contribution on (Hi + H 2 ) 2 (see, for instance, Ref. Q). Expanding in (Hi ± H 2 ) 2 and taking the 
limit Hi — > H 2 = H, we arrive at the desired result. 

All of the above results can be also obtained in a purely thermodynamical framework using the small corrections 
theorem [[l3|. Indeed, the latter states that small corrections to all thermodynamic potentials are the same provided 
that the appropriate thermodynamic variables are kept constant. Therefore, we write, 



s h f = s H n (p (H) , H) 



(A10) 



an no variation of p (H) is done in the r.h.s. This is in correspondence with eq. (A9) where only the level density is 
varied with the field. Expanding around (p), we find 



<hi F ^ SniU.il'} . U) -8 B (—«/,), H) ) (// - (//}) 



which, again, is in complete correspondence with the expansion in eq. (A6). Noticing that 

N(p(H),H)=N 



(All) 



(A12) 



we reproduce eq. (A9) (neglecting the difference between the thermodynamic and si ngle -particle densities of states, 
dN/dp and (p) respectively, and after some manipulations analogous to those in eq. (A7)) as follows: 



SF H = S H F-S H n((p),H) 

w -5 B (N((p) ,H)-N(p (H) , H)) (p - (p)) 
„ dN((p),H) 



dp 



(S H (p-(p)))(p-(p)) 
1 



(A13) 



S H (N - (N)Y 



APPENDIX B: PERSISTENT CURRENT IN THE PERTURB ATIVE APPROXIMATION 



Using cqs. (0), (JT9j) and ([14]) we find, upon reducing the integral to the Matsubara sum, 



Ipert (<p) 



t 2 d 



E 



2 ^ 2{p) 9^ 1 J mT+(27rrHr a 
With the help of the identity u~ 2 = J °° t exp (— ut) dt the latter can be re- written as 



(Bl) 



Ipert (4>) 



E 



1 dr H ' - 



2/3tt 2 2T (p) 



TlX 



n \ sinli ( TT.c ) I ^ C T 



(B2) 
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which coincides with the respective equation^] in Ref. Q . From eq. @ it follows that § 

1 dr^ 1 r / x \2 , s 2 8<?(2)e£: c <j 



L pert 



Lpert 



(0) 



2/3tt 2 2ttT (p) 6*^6 7 



<6 L VsinhaJ X 2Bn 



28tt 2 T(p) O : 



1 dr^ 



2/3tt 2 2T (p) 



2/3tt 2 2tt (p) 



< < 00 



where ^ (2) = 7r 2 /6 and C is found from the condition t^ 1 ~ T and is given by 



T 



4ir 2 E r 2 



< 00 



(B3) 
(B4) 

(B5) 



For completeness, we also derive the result accounting for the Aharonov-Bohm periodicity in narrow rings. This is 
easily accomplished by using the following expression for the perturbative level-density correlation function 



R pert (lu,t h x ) = Re 



1 



An 2 E r 



2(1) 



(B6) 



As a result, we find the following expression for the persistent current 



'-pert 



{4>) 



t 2 d 



oo oo 

2 in 1 lip) ,)<> E E 

vii n— — oc m—1 



mT + 2ttE 



(B7) 



which can be evaluated using and identity 



E 



7r sinh (27r6) 



\b 2 + (n + a) 2 cosh (27r6) — cos (2-7ra) 



(B8) 



(notice that for small a and b the s um i s dom inat ed by the n = term - in agreement with the procedure adopted in 
the bulk of the article). From eqs. (B7) and (B8), we find 



Ipert 



(</>) 



vrT 2 d d 



2/3?r 2 2 (p) E c dT , 



E 



i 



sinh 



(a/5 



,, , J2g m cosh (y^m)- cos (^) 



(B9) 



In Ref. jjj, the sum appearing in this equation was evaluated using the Euler-Maclaurin method. However, the latter 
does not work well for small values of the argument. Consequently, here we choose a different approach. Namely, 
we use the fact that 2itT <gC E c and expand the argument for m <C E c /2nT and extend the sum to infinity thus 
neglecting the correction of the order 2irT/E c 



Ipert (0) 



} _, _ T sin (4^0/0 o ) (z) + \z^ [z)\ 



E, 



z = — [1 - cos (4tt0/0 o )] 

TTl 



(BIO) 
(Bll) 



where tp^ ( z ) is the polygamma functio n. U sing ip^ (z) + \zip ^ (z) $ (2) and ip^ (z) + ^zip^ (z) z ^^° (2z) 1 , 
we recover - in the limit <C 0o - eqs. (B3) and (B4) from eq. (BIO). 



2 Notice that, due to a misprint, the factor T 1 is missing in eq. (18) of Ref. Q. 
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